Inverse problem and Bertrand's theorem 
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The Bertrand's theorem can be formulated as the solution of an inverse problem for a classical 
unidimensional motion. We show that the solutions of these problems, if restricted to a given class, 
can be obtained by solving a numerical equation. This permit a particulary compact and elegant 
proof of Bertrand's theorem. 
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£^ ' In classical point mechanics, since central potentials are isotropic they're at the basis of every two body interactions. 
. Among them we find in particular Hooke's and Newton's potentials which possess very specific properties. A part of 
them were already known by these pionners of modern theoretical physics. For example, they're the only potentials 
t— I . presenting elliptic bound states'orbits and the ellipses associated to each of them are in dual relation0. 

The question to determine if others central potentials could generate closed bound states orbits for every values 
of the initial parameters of the motion (energy and angular momentum) stayed opened for almost two centuries. In 
1873, J. Bertrand has shown that the answer to this question is negative. Since thirty years, a number of others proofs 
^vq . of this fundamental result have been proposed. With an exception 0] , [3] , every proof schemes can be decomposed in 
7— I ' three steps. The first one establish that for the researched potentials, the angular period, also called apsidal angle, is 
necessarily independent of energy and angular momentum. In the second step, one shows that constancy with respect 
to angular momentum, when applied to orbits closed to circular ones, leads to retain only power law potentials. Both 
of these steps are common to every proofs. Only the last step is handled differently. It consists to show that on! 
Qh' Hooke's and Newton's potentials lead always to closed orbits. Most of the proofs use a perturbative approach 
r- 1 • ,11,01,0, 0- The original proof , [l(| and the inspired ones [ll[ follow rather a global approach. Since the problem 
is to determine a class of potentials from informations about the orbital angular period, the most natural and direct 
approach is to interpret it as an inverse problem. This approach has already been used by Y. Tikochinsky |8| and in 
a different and more general setting by E. Onofri and M. Pauri [lij. In Tikochinsky's work 0, it effectively runs into 
a simplified proof, but the inverse problem treatment still necessitates tedious perturbative calculations. 

In this paper, we propose to treat the inverse problem from a new point of view, then obtaining a particularly 
compact and elegant proof of Bertrand's theorem. Moreover, we will show that this formulation permits as well to 



' deal with the problem perturbatively. 

Basics about the motion of a particle in a central force field 

Let's consider a particle submitted to a central force field deriving from the potential U (r). The angular momentum 
L (t) = |?o x t a is then conserved and every orbit lies in the plane perpendicular to L . 
In polar coordinates (r, ip), the equation of motion in this plane gives 11 1 , , [H3] , 14\ : 



r (t) + £V£(r(t)) = 

v (*) = ^kt) 



(1) 



where Vl(t) is the radial effective potential, that is the sum of the initial potential and the centrifugal barrier 

L 2 
2mr' 2 



T 

term 2 3 which intensity depends of the angular momentum 



V *)= U W+-£* ( 2 ) 
The radial coordinate r(t) describes an autonomous unidimensionnal motion, those of a particle submitted to Vl{t). 
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V L (r) 




The solution r(t) of equation |T]) is given, at least implicitly, by Barrow's formula : 

^ dp 



t 



(?)' 



VE-Vl(p) 



(3) 



where E is the conserved energy of the system and where we have chosen the initial conditions to — 0, r(0) = 
r , ip(0) = ipo- 

The angular coordinates (p (t) is obtained from r(t) by a simple integration : 



V (*) - = — / 



eft 



(4) 



We therefore lead to a complete parametric description (r(t), ip(t)) of the motion with respect to the time. 
Let's introduce some elements of the vocabulary usually used for the description of this type of motion. 



= infr^ such as E = Vl(t"j) is called a pericentral radius and we have E < VL(r m — e), V e > sufficiently 



small. Every point A m such that 



OA n 



= r m , is a pericenter. 
a apocentral radiu 

tm is an apocenter. 



* r M = supr, such as E — Vl(^) is called a apocentral radius and we have E < Vl{t m + e), V e > sufficiently 

i 

small. Every point A M such that OA M 
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FIG. 2: Radial effective potentials for power law potentials 



* r m and tm are the apsidal distances of the motion the vector positions ~r A m & ~~t a m °f pericenters and 
apocenters are called apsidal vectors. The angle $ between two consecutive apsidal vectors is the apsidal angle. 

If ru < +00 , the orbit is bounded. If therefore r m > 0, the radial motion is an oscillatory motion between the 
two extremal values r m and tm- The orbit is then localized in the annulus r m < r < tm ■ Such a bounded orbit will 
only be closed at the condition that 2$ (which is the angle between two consecutive pericenters or two consecutive 
apocenters) be commensurable with 27r, that is : 

$ = -7T e ttQ 

q 

where p and q are integers. In this case indeed, after p revolution around the origin, the particle makes q radial 
oscillations. 

In every others cases, that is $ ^ Qtt, the bounded orbit is everywhere dense in the annulus r m < r < tm- This 
kind of orbit is called a rosette. 



Circular orbits 

If the orbit is bounded there exists at least one absolute minimum for Vl{t) on the interval [r m ,rM]- We'll note R 
the corresponding value of the radial abscises and : 

V R = V L {R) = U{R) + ^ (5) 

The minimum conditions give : 

/ Vl(R)=0^^=R 3 U'(R) 

\ V£(R) = ^"(R)+3u'(R) > Q W 

If E = Vr, the only authorized value for r is R and : r(t) = R = r m = r M , Vi with r(<p) = R, Vtp, which corresponds 
to a circular orbit with radius R. 
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Along such a circular orbit we have : 



r(t) = R 



(7) 



If the initial potential U(r) is such that the function f(r) — r 3 U'(r), is locally bijective around each minimum of 
Vl(t), the first of the above conditions © shows that we can indifferently choose to characterize circular orbits by 
their radius R or by the associated angular momentum L. 

Let's finally note that in order to admit a finite circular orbit, the power law potential U(r) ~ r v must satisfy v > 
or — 2 < v < 0, that is U(r) must be of the form : 



U v {r)=kr\ u>0 
U- V {r) = -kr- v , 0<v<2 

In each case we have a unique circular orbit with a radius respectively equal to 



(8) 



u V vkra j 

R - (J±- 
IX - V ~ \ vkrn 



V > 

< v < 2 



(9) 



Clairaut's variable and Binet's equation 



If our only ambition is to describe the orbit, it is possible to obtain its polar equation in a direct way. To do it, it's 
sufficient to note that, if we consider now r as a function of <p, eequation (JTJ becomes (j^ = m Jl^ -^j '■ 

-d^-^{—) + — c3-^(r(v)) = (10) 

Its solutions furnish directly the orbital equation under the form r = r(<p). This equation is considerably simplified if 
we do the following change of variable x = (due to Clairaut [3|)- Clairaut's variable x, is nothing else, upon 
to a constant factor, that the inverse radial variable, x and tp are called Clairaut's coordinates. 
With this change of coordinate, equation (fit)]) then becomes : 

' ~W' L (*(*,)) = (11) 



dtp 2 m 



where : 



W L (x) = V L (— ) =Knx 2 + u(— ) (12) 
\mx ) 2 \mx ) 

This identity is more known as Binet's formula or Binet-Clairaut's equation. 

When tp varies, Clairaut's variable x then describes a one-dimensional motion those of an effective particle of masse 
m submitted to the Binet-Clairaut's potential Wl(x). The evolution parameter of this motion is the angular 
position (p, growing with time. When r(tp) makes an oscillation between the values r m and rjj then x((p) makes a 
corresponding oscillation between the two extrema a;> = L and x< — L . Let's note that to each extremum of 

Vl{t) in r = R corresponds a extremum of the same type for Wl(x) in xq = Concerning the curvature of the 
Clairaut's potential near a minimum xq, that is a circular orbit, it writes : 

RU"(R) + 3U'(R) ^ n 
W (x ) = m > (13) 

Choosing the angles origin at an apocentral vector ~?a m and taking the initial condition x (0) = x < , Barrow's 
formula, when applied to Binet-Clairaut's equation, gives an implicit solution for the orbital equation : 



(14) 
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Then the apsidal angle <I>, which is Clairaut's motion and the half-period of the radial oscillation is : 

fm\h f x> dx , . 

^ £ )-(t) L T^wm (15) 

This expression takes a very compact form if we use the semi-derivative's concept [2fJ]. The semi-derivative is the 
integral operator defined as [15| : 

i r E , i 



Dlg(E) = -= / dw g'(w) (16) 

V 77 Jv R v» — w 

where the function g(E) satisfies g{Vn) = 0. 

The semi-derivative is a particular case of fractional derivative, notion about which an abundant mathematical 
literature is available S3 , S3 , S3 and which has found today many physical applications S3 , S3 , [2(| ■ 

Putting Ax(w) = x> (w) — x< (w), where x> (w) and x < (w) are the reciprocals of Wl(x) defined on each branches 
of this last, on both sides of xo, equation (fl"5|) becomes : 

/ TV) 7T 1 

$(E,L) = ^J— D»Ax(E) (17) 

If the apsidal angle's notion becomes meaningless for a strictly circular orbit, we can however try to calculate the 
limit value it takes when the considered orbit is in the neighborhood of a circular one. From the point of view of the 
radial or Clairaut's motions (r(t) or x(ip)), this corresponds to small oscillations near the equilibrium values R and 
x = of potentials Vl{t) and Wl(x) respectively. 

Let's recall [13, S3 

that in the small oscillations limit, a particle whose potential W(x) has a non zero curvature 
in the vicinity of the equilibrium position xo is, at the first order of approximation an harmonic, (that is isochronous) 
oscillator with frequency : 



<W"{x ) 

uj = y 

In the present context, this simply gives (fT5|) 



(18) 



where the angular period of radial oscillations in the vicinity of a circular orbit <I>c(i?), is now independent of the 
energy E. 

Knowing the function <f>c(i?), the above identity becomes a second order linear differential equation for U(R) whose 
solution is readily obtained as : 



„, -r'^(3-(w) 2 ) 



U(r) = J dr'e P V voiwy } ( 20 ) 

In the following we will specially refer to the case where $c(R) is a constant <E>c, independent of the radius R (that 
is of the angular momentum L of the particle). In this case, from the above formula we obtain two distinct functional 
forms for U(r) : * If <I>c ^ defining the exponent v as $c = 1 v ^ ^ ' 

U(r) = Ar v + B (21) 

* If $ c = (that is for v = 0) : 

U(r) = Ahir + B (22) 

Therefore, there are only two types of potentials for which the apsidal angle $c of a bounded orbit near a given 
circular one, is a constant, independent of the characteristic parameter L : power law potentials and logarithmic 
potentials. 
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Inverse problem for the Clairaut's potential 



The inverse problem for a classical one dimensional oscillator consists in determining the oscillator potential from 
the variation law of the period as a function of the total energy. Then to determine Clairaut's potential Wl(x) from 
the variations of the apsidal angle &(E, L) as a function of E is a problem of this type. Fractional integro-differential 
calculus as introduced above gives a very direct way to the solution [20| . Indeed, the semi-derivative operator Di 
(fT6|) admits, on the set of bounded functions near Vr, an inverse called semi-integral 15|,[2(j and defined as : 

D~*g(E) = A= T dw J. g(w) (23) 

V n Jv R v E -w 

For our purpose, we'll more particularly note the following result for the semi-integral of the constant 1 : 

D^l = -y/w-VR (24) 
7T 

Since Ax(w) satisfies the required condition, equation (|17|) linking Ax to <f>(E,L) is immediately inverted and the 
inverse problem's solution for Clairaut's motion is simply given by the implicit formula (2p| | : 

Aa; (E) = x>(E) - xAE) = J —DZ?$(E,L) (25) 

V mix 

Unfortunately, this equation yields only the difference x>(E) — x<(E) for a given &(E,L), which doesn't generally 
allow to determine uniquely the two branches x > (w) and x < (w) of the potential's reciprocal function. For a given 
Ax (w) there exists an infinity of possible multiform x (w). Therefore, if our goal is to precise what potential function 
leading to a given $(£', L), we see that this problem will generally admit an infinity of solutions [13j |. 

To determine Wl(x) uniquely it's necessary to add supplementary constraints, that is to limit the research to a 
more restricted class of potentials. For instance, if we consider only potentials which are symmetrical with respect to 
the axis x — Xq, we obtain : 

1 _i 

x> (w) = D w 2 $(w,L) +x = 2x -x<(w) (26) 

V 2r7J7T 

But even if we have obtained x > (w) and x < (w), we still have to invert them to find an explicit form of Wl{x). 
Usually this step is not analytically tractable. 

We can avoid this difficulty by translating the identity ([23]) into a functional equation for Wl{x) [2(|. Indeed we 
have : 

W L {x) =W L (x + Ax) (27) 

in every point x of the interval I — ]0, Xq], when Aa; (w) is calculated in w = Wl(x). 
Using ([23)) . we then obtain the following functional equation for Wl(x) : 

W L (x) = W L \x+J—D^$(W L {x),L)) (28) 

In the special case where $>(E,L) — $c(E) constant with respect to E ("isochronous motion"), this gives (|24p : 

W L {x) = W L (x + a L ^W L (x)-V R ^ , Vxel (29) 

with a L = X /I**2i&. 

A priori this equation doesn't seem particularly simple to solve analytically. However if we suppose we can limit our 
research to a restricted class of potentials W(cc; {vi}) with a given functional form but depending upon m parameters 
i v i)i—i m i we are leading to a system of n numerical equations for the fj, with n arbitrarily large : 



' W{x 1 -{v 1 })=w(x 1 + ^dJ<S>(W(x 1 -{v 1 }),L)-{u 1 }) 
W{x n - { Vi }) =w(x n + ^D-^(W(x n ;{^}), L); {^}) 



(30) 
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where the Xj are n points on /. 

Specially, for a family of isochronous potentials depending upon a unique parameter v [m — 1), we obtain the 
following equation : 

W{x x ;v) = W (x x + a L y/WUx^v)- V R ; i/) (31) 

where x\ is arbitrarily chosen on / (n = 1). 

We can equally chose for x x the limit value + . With regard to the possible behaviors of Wl(x) at this point we'll 
consider only the two following cases : 

* Wl(x) tends to a finite value Wl,o = Wl(0 + ) in (this corresponds to the case where U(r) tends to a finite value 
when r — > oo). Equation (|31[) becomes : 

W Lfl (u) = W L (a L tJWlvW-Vr; uj (32) 

* Wl(x) diverges as Cox~^° when x ^ and as CqoX^ 00 when x — > oo (/io, Moo, Co, > being a priori t> 
dependent). This corresponds to the case where £/(r) diverges as r^ when r — > oo and as r - ^ 00 when r — * 0. 
Equation leads then to the relation : 

CaT" ~ OoC^a^x-^ ( ^ ' = 2 2 (33) 
The Bertrand's theorem 

All the preceding results permit to obtain a particularly compact proof of Bertrand's theorem. This last establish 
that : 

* The only potentials for which every orbits near a circular one are closed independently of energy 
and angular momentum are the Newton potential and the Hooke potential. 

If one except the proofs based on the necessary existence of supplementary constants of motion, all the other 

proof's schemes can be decomposed in three steps, the two firsts being common to all proofs. 

First step : Obtaining the variation law for the angular period $ as a function of E and L 

The first step is the keystone of the proof. In order that all the bounded orbits be closed it's necessary that on 
all the intervals of admissible values for the parameters L and E the corresponding values of the apsidal angle $ are 
rational : 

$(E,L) = -7T e ttQ 
q 

that is belong to a discrete set. 

$(-E, L) being supposed to vary continuously with L and E, the fact that its image is contained in a discrete set 
implies that it's a constant with respect to L and E. 

If we note $c the limit apsidal angle apsidal near the considered circular orbit, then we necessarily have : 

&(E,L) = $ c SttQ (34) 

for all the considered values of E and L. 

This extremely constraining constancy condition is the source of the result. 

Second step : Selection of the potentials for which the limit apsidal angle near a circular orbit is not L dependent 

The second step consists to use the preceding results (|2"Tj) . (|2"2"j) concerning the potentials possessing circular orbits 
near which the limit apsidal angle is L independent. Among them, the logarithmic potential Ui og (r) = fclnr + 
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cste, k > 0, can be immediately excluded since it leads to a limit apsidal angle <&c = 3^ which is not a rational 
multiple of it and then doesn't satisfy the condition (|34)) . 

Therefore, the only permitted potentials U(r) are of the type U± w {r) — ±kr ±l/ , k > 0, and their associated 
Clairaut's potentials writes (O : 

W ±v {x) = \m{x 2 ±A ± x^), A ± = 2-(-) (35) 
2 m \m J 

with v > for the + sign and < v < 2 for the — sign. 

In both cases, one has a unique circular orbit for which the limit apsidal angle is (|21[) : 



Third step : Determination of the admissible power law potentials 

Now it remains finally to select among all the possible values of the exponent v, those permitting to satisfy the 
required conditions. The differences between the various proof's scheme appear at this level. The great majority lie 
on a local study in the circular orbit's neighborhood, the corrections be ing determined by an adapted perturbative 
approach The original proof [§], [13] and Arnold's one proceed from a global approach. Now, 

we are going to see that our inverse problem's formulation as developed in the preceding paragraphs permits to solve 
the asked problem in a very simple and direct way if we consider globally the potential behavior. We will show that 
this formulation permits equally a local perturbative treatment, requiring however tedious calculations, which is a 
common feature of this kind of approach. 

Let's come back to the asked problem. We have to determine all the (Clairaut's) potentials leading to a given 
variation law for the (angular) period as a function of energy. As noted by Tikochinsky, we are typically in the frame 
of an inverse problem [8j]. Since in our case, the variation law is a constant, we are therefore leaded to determine 
among all the Clairaut's potentials in the families (W ll (x)) l/>0 et (W- v (x)) 0<l/<2 (|3"5|) those being isochronous, that 
is satisfying the functional equation (|2U|) . 

Global approach As seen previously, by choosing an adapted value for x%, the functional equation (|29[) results 
in numerical equation for v (|3"Tj) . Let's take for x\ the limit value + . The behaviors at the origin of the potentials 
W v {x) and W- V {x) are : 

WJx) ~ \mAx~ v -> +oo 

> > 

W- V {x) -> W L>0 = 



i40 



* For (W- u (x)) Q<iy<2 , equation gives : 



W. v (aVWnl) = k ( ±- V" 2 f 2 - f ~) " 2 ) = (38) 



v vkm I V V-4 
where in the present case a = \ — '^ £ - = \ — , 2 . 

* y m 7r y m \J2±v 

We then obtain a transcendental equation for v : 

e -5 ln ^=2 (39) 

whose solutions are readily obtained as v = 1 and v = 2. 

Since < v < 2, only v — 1 is an admissible value. This corresponds to an initial potential U(r) = —-, that is the 



Newton's potential. 

* For (W„(x)) ly>0 , identity ([33"]) takes the form : 



1 = (40) 



with Coo = \m and a = 2 J mi % +v) ■ Then 



= 1 v = 2 

2 + v 



W L (x) 



w 



Vr 




£+ 



FIG. 3: Lateral displacements 



In this case, the initial potential U(r) is harmonic : U(r) = kr 2 . 

To summarize, the only potentials satisfying the required conditions of the Bertrand's theorem are the Newton 
potential and the Hooke potential, which achieves our proof of the theorem. 

Perturbative approach As mentioned previously, it's possible to recover this result via a perturbative resolution 
of the functional equation (|29|) satisfied by Clairaut's potential. 

For that purpose, let's introduce the left and right lateral displacements e_ and e+. They measure the respective 
distances between Xq and the two branches of the potential function Wl (%) ■ 



E- = Xo — X 

e+ = Ax — E- 



(41) 



Then, the functional equation satisfied by Wl(x) is in every point x < xq : 

W L {x - e_ ) - V R = Wl (x + e+) - V R (42) 
If we consider an isochronous Clairaut's potential Wl(x), the relation between the lateral displacements writes : 

e+ = Ax - s- = a^Jw L (x - eJ) - V L , m - £- (43) 

where a = \ I — . 



For small amplitudes displacements, we can expand the left side of the above identity in power of e_ 
First, we have : 



W L {xo-£-)-V R 



ImcV (f; (-l) n a n eA 
\n=0 / 



(44) 
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where W^(xq) = mcu 2 and : 



1 2Wi n+2 \x Q ) 



(n + 2)! muj 2 

Inserting this expansion in the relation (|43|) between lateral displacements, we obtain : 



(45) 



1 - a l£ - + a 2 e 2 _ + ... + (-1)™ a n e n _ + ... - e_ (46) 



- 1 



(2n - 3)!! +1 / , p x \™ 

— E 2 n ra! £ + + - + (-1) P a P e P - + ■■■) (47) 



7T 

2^of (->"- ■> '!!.„ + , 

n=l 

that is, up to the fourth order in e : 

e+ = (2 7 - 1) - + e 3 _ 7 (oa - f ) + ei 7 (-a 3 + ^ - f ) + O (e 5 _) (48) 

with 7 = 

Then, if we expand both members of the functional equation (|42[) in power series of e_ and e+ respectively, we hnd 

e 2 _ (1 - + a 2 £* + ... + (-1)" a n e n _ + ...) = e\ (l + ai£+ + a 2 £+ + •■■ + a n e+ + ■■■) (49) 
Inserting in this identity the above expression of e+ 1481 we arrive at the formula : 

1 - a l£ - + a 2 e 2 _ = (2 7 - l) 2 + (2 7 - 1) ((2 7 - l) 2 - 2 7 ) ai£_ 

+£ 2 _ (l 2 a\ + 2 7 (2 7 - 1) (a 2 - ^ - 3 (2 7 - l) 2 7 a 2 + a 2 (2 7 - I) 4 

Identifying order by order the coefficients in each side we then obtain : 

7= fiiSa = 1 

0.ai = (50) 
a 2 = \a\ 

The first of these equations simply translates the fact that, in the small oscillations limit, x(ip) executes harmonic 
oscillations, which period is : 

2tt , s 

2$ c - . 51 



The second equation tells us any information. As to the last, it connects the third and fourth derivatives of the 
Clairaut's potential in x : 



r(4) . . 5 



(wf(x ) 



W L M = o ( 52 ) 

Of course, these constraints on the Clairaut's potential, generated by the isochronism condition, are exactly the same 
that those obtained in a singular perturbation expansion [7| , due to the suppression of the corrections on frequency. 
Therefore to generate isochronous oscillations near Xq (period 2&c — , 2?r ), the Clairaut's potential has to 



satisfy this constraint. 

If we refer to the class of potentials W±„(x) ([35]), this gives [W$(x ) = (=F ^~ 3) wjfj (^o), W£ v {x Q ) = m + 
v{v±l) A±x^ u - 2 = m (2 ± v)) : 

y xq 3 ui (2 i is j I 
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For a potential W u (x), v > 0, we obtain : 



2{2-v) 
3x a 







(54) 



which implies v = 2. 

For a potential ^-^(x), < v < 2, we have on the other hand 








(55) 



which implies W_„(xq) = that is u = 1. 

We recover the preceding result. Nevertheless, it has required much more tedious calculations, which, as we already 
pointed out, is an inherent feature of the perturbative schemes of proof [1] , [H , 0] , Q • 
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